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ABSTRACT 

We have analyzed the far-nltraviolet (FUV) spectra of six elliptical and 
SO galaxies in order to characterize the stellar population that produces the 
ultraviolet flux in these galaxies. The spectra were obtained using the Hopkins 
Ultraviolet Telescope (HUT) during the Astro-2 mission aboard the Space 
Shuttle Endeavour in March 1995, and cover the spectral range from 820 to 
1840 A with a resolution of 3 A. These data, together with the spectra of 
two galaxies observed with HUT on the Astro-1 mission, represent the only 
FUV spectra of early type galaxies that extend to the Lyman limit at 912 A 
and therefore include the “turnover” in the spectral energy distribution below 
Lyman alpha. 

Using an extensive new grid of LTE and non-LTE synthetic spectra which 
match the HUT resolution and cover the relevant parameter space of temperature 
and gravity, we have constructed synthetic spectral energy distributions by 
integrating over various predicted stellar evolutionary tracks for horizontal 
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branch stars and their progeny. When the compnted models are compared 
with the HUT data, we hnd that models with snpersolar metal abnndances 
and helinm best reprodnce the flnx across the entire HUT wavelength range, 
while those with snbsolar Z &: V fit less well, partly because of a signihcant 
flux dehcit shortward of 970 A in the models. High Z models are preferred 
because the contribution from the later, hotter, post-HB evolutionary stages 
makes up a higher fraction of the sub-Lyman a flux in these tracks. We hnd 
that AGB-Manque evolution is required in all of the hts to the HUT spectra, 
suggesting that all of the galaxies have some subdwarf B star population. At 
any Z & U, the model spectra that best match the HUT hux are dominated 
by stars evolving from a narrow range of envelope mass on the blue end of the 
horizontal branch. 

The Astro-1 and Astro-2 data are also the hrst with the resolution and 
signal-to-noise needed to detect and measure absorption lines in the FUV 
spectra of elliptical galaxies, allowing a direct estimate of the abundances in the 
atmospheres of the stars that produce the UV hux. We hnd that most absorption 
features in the spectra are consistent with Z = 0.1 Zq, signihcantly lower 
than the abundances implied by the best-htting spectral energy distributions. 
However, given the strong observational and theoretical evidence for dihusion 
processes in the atmospheres of evolved stars, the observed atmospheric 
abundances may not rehect the interior abundances in the population producing 
the ultraviolet hux in elliptical galaxies. 


Subject headings: galaxies: evolution — galaxies: abundances — galaxies: 
stellar content — ultraviolet: galaxies — ultraviolet: stars 


1. INTRODUCTION 

An understanding of the stellar populations in elliptical galaxies will shed light upon 
conhicting theories in cosmology, galactic evolution, and stellar evolution. However, 25 
years after the initiation of space-based observations and the consequent extension of data 
into the ultraviolet (UV), there are still considerable uncertainties regarding the chemical 
composition in elliptical galaxies and the evolution of their stellar populations. 

The spectra of elliptical galaxies and spiral galaxy bulges exhibit a strong upturn 
shortward of 2700 A, dubbed the “UV upturn.” Characterized by the (misso — V) color, 
the UV upturn shows strong variation (ranging from 2.05-4.50 mag) in nearby quiescent 
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early-type galaxies (Burstein et al. 1988). The strength of the UV upturn is positively 
correlated with the strength of Mg 2 line absorption in the V band, in the sense that the 
(^^1550 ~ V) color is bluer at higher line strengths, opposite to the behavior of optical 
color indices (Burstein et ah 1988). Opposing theories have been devised to explain this 
correlation. In one camp, Lee (1994) and Park & Lee (1997) suggest that the UV flux 
originates in the low metallicity tail of an evolved stellar population with a wide metallicity 
distribution. In the other camp, several groups (Bressan, Chiosi, & Fagotto 1994; Greggio & 
Renzini 1990; Horch, Demarque, & Pinsonneault 1992) propose that metal-rich horizontal 
branch (HB) stars and their progeny are responsible for the UV flux. These different 
metallicity scenarios in turn argue for different ages for the stellar populations in these 
galaxies. Ages exceeding those of Galactic globular clusters are required in the Park & Lee 
(1997) model, while ages as low as 8 Gyr are allowed in the Bressan et al. (1994) model. 

In these two scenarios, the EHB stars are drawn from either tail of the metallicity 
distribution. However, it is also possible, indeed perhaps more likely, that the EHB stars 
arise from progenitors near the peak of the metallicity distribution (cf. Dorman, O’Gonnell, 
& Rood 1995), but represent a relatively rare occurrence. Observations of a bimodal HB 
within the metal-rich open cluster NGG 6791 suggest that this might be the case (Liebert, 
Saffer, & Green 1994). If the same mechanism is at work in both the elliptical galaxies and 
in NGG 6791, the correlation of (misso — V) with the global metallicity of the galaxy might 
indicate that this rare path of stellar evolution becomes less so at high metallicities. 


1.1. Horizontal Branch Morphology 

How can two diametrically opposed interpretations arise from the same observational 
data? The answer lies in the parameters that govern HB morphology, such as metallicity. 
In general, for stars of solar metallicity or less ([Fe/H] < 0), higher metallicities yield redder 
HBs and lower metallicities yield bluer HBs, i.e. bluer HBs have more stars on the blue 
side of the RR Lyrae instability strip, and redder HBs have more stars on the red side. 
Metallicity is the implied “first parameter” of HB morphology. In his spectrophotometry of 
13 Milky Way globular clusters, Gregg (1994) gives an excellent observational example of 
the classical effect of metallicity on HB morphology; Lee, Demarque, & Zinn (1994) provide 
theoretical examples. However, it has been known for decades that other parameters affect 
HB morphology (cf. Lee et al. 1994 and references therein). In this way, the metallicity 
debate in elliptical galaxies is tied to the “second parameter” debate that has been raging 
in the study of globular clusters. 

Helium abundance also plays an important role in the stellar structure of zero-age HB 



(ZAHB) stars (Dorman 1992). These stars are characterized by a helium-rich core and a 
hydrogen-rich envelope. At the core/envelope boundary, the density p and mean molecular 
weight /i change such that p/p, is continuous. As the helium abundance (D) in the envelope 
increases, the p gradient across the core/envelope boundary decreases, which in turn 
decreases the temperature gradient across this boundary. Shallower temperature gradients 
produce higher shell temperatures, and the nuclear reactions in the hydrogen-burning shell 
are highly sensitive to the temperature. At the red end of the ZAHB, the shell luminosity 
in these massive envelopes makes a signihcant contribution to the total luminosity, and 
so an increase in temperature results in a signihcant increase in the total luminosity, 
but does not change the envelope structure signihcantly. At the blue end of the ZAHB, 
the shell luminosity in these less massive envelopes makes a smaller contribution to the 
total luminosity, and so the total luminosity does not change signihcantly with increasing 
temperature, but the envelope structure does. Thus, an increase in Y on the blue end of the 
ZAHB signihcantly increases Tefr and slightly increases the luminosity, while an increase in 
Y on the red end of the ZAHB signihcantly increases the luminosity and slightly increases 
Teff. Across the entire ZAHB, an increase in Y increases the rate of nuclear reactions, 
resulting in more rapid consumption of the envelope. 

Lee et al. (1994) argue that age is the second parameter driving HB morphology (see 
also Sarajedini, Lee, & Lee 1995). They also investigate the role of He abundance, CNO/Fe 
abundance ratio, and core rotation, but rule out these latter candidates for the second 
parameter because of conhicts with observations. Lee et al. (1994) hnd that age variations 
of several Gyr can produce dramatic diherences in HB morphology, in the sense that young 
HBs (age = 11-12 Gyr) tend to be red, while older HBs (age = 15-16 Gyr) tend to be 
blue. This is because the mass of a red giant at the helium hash decreases as the age of a 
stellar population increases, which in turn decreases the mean mass of the HB stars (Lee et 
al. 1994 and references therein). Although this scenario is true in any given system, it must 
be stressed that when comparing diherent populations, the total mass loss on the red giant 
branch (RGB) may not have the same mean value and dispersion, and thus diherences in 
HB morphology cannot necessarily be attributed to diherences in age alone. When more 
than one parameter drives the HB mass distribution, it may be difficult to separate age 
from these other ehects. 

For an HB with a single metallicity, the envelope mass {Menv) determines the color 
location of stars on the HB. According to theory, stars with lower envelope mass are bluer 
(hotter) and stars with higher envelope mass are redder (cooler). Stars evolving oh the 
HB at core helium exhaustion will then follow one of three diherent evolutionary scenarios, 
again depending upon the envelope mass. Dorman, Rood, & O’Gonnell (1993, hereafter 
DR093) explain these evolutionary scenarios in detail; we will briehy summarize them here. 
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The first scenario occurs for an HB star with M^nv ^ 0.07 Mq (in the case of solar Z 
and Y). Such a star evolves along the asymptotic giant branch (AGB), and after reaching 
the thermally pulsing (TP) stage, much of the envelope is removed and the star no longer 
sustains a deep convective exterior. At this point it evolves rapidly to higher temperatures 
at hxed luminosity, and then fades along the white dwarf (WD) cooling track. An example 
of a classical post-asymptotic giant branch (PAGB) evolutionary track (Vassiliadis & Wood 
1994) is given in Fig. 1, along with the integrated far-ultraviolet (FUV) spectrum of a 
population of stars entering that track at the rate of one star per year. An AGB path 
that might precede the PAGB evolution is also shown (DR093), but this cool, red portion 
of the evolution (from HB to TP) makes a negligible contribution to the integrated FUV 
luminosity in the HUT wavelength range. 

The second scenario occurs if an HB star has enough envelope mass to reach the 
AGB, but not enough to maintain a convective exterior all the way to the TP stage 
(0.03 Mq Menv ^ 0.07 Mq for Y = Yq and Z = Zq). These stars, dubbed post-early 
AGB (PEAGB; Brocato et ah 1990), evolve away from the AGB before reaching the TP 
stage, moving to higher temperatures at a constant luminosity which is less than that in the 
PAGB scenario. Because of the lower luminosity, and because models with lower envelope 
masses leave the AGB with more remaining hydrogen fuel, the PEAGB stars have longer 
lifetimes. Thus they emit more UV light over those lifetimes, as compared to the stars 
following the PAGB tracks. Fig. 1 shows an example of PEAGB evolution (DR093) and the 
corresponding FUV luminosity. This particular PEAGB track exhibits helium shell flash 
TP “loops” which take place on the AGB if the envelope is more than a few hundredths of 
a solar mass. Although they are dramatic in the HR diagram, these short-lived features do 
not signihcantly alter the integrated spectrum. 

The third scenario occurs for those HB stars with such small envelope masses 
{Menv ^ 0.03 Mq for Y = Yq and Z = Zq) that they are unable ever to develop deep 
outer convection zones, and thus are unable to reach the AGB at all. In the shell burning 
phase they are referred to as AGB-Manque (AGBM) stars, and they spend 10^ yr at 
relatively high luminosity (~ 100 Lq) and high temperature (Tefj > 30,000 K), after which 
they evolve directly to the WD cooling track. The amount of FUV light emitted from 
AGBM stars is much higher than that emitted from PAGB or PEAGB stars, as shown 
in Fig. 1, again due to their longer evolutionary time scales. Following the nomenclature 
of DR093, we will use the term extreme HB (EHB) for those HB stars that evolve along 
PEAGB or AGBM evolutionary paths. 

The different evolutionary paths a star may follow after leaving the HB are affected by 
the metallicity and helium abundance of the star. By itself, the envelope opacity at high Z 
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has a tendency to “damp down” the hydrogen burning shell, so that high-metallicity stars 
are fainter (Dorman 1992). It is usually supposed that enhanced metal abundance is seen in 
stellar populations accompanied with increases in the helium fraction, according to a simple 
linear formula: AY/AZ = constant. Values for this constant are, however, determined 
from observations close to Z = 0 rather than Z = Zq, since what is often sought is the 
primordial helium abundance for cosmological studies. There is little information about 
the helium abundance of metal-rich populations and the subject is quite controversial (see 
Dorman et al. 1995 and references therein; for a recent study, see Minniti 1995). For high 
He envelope abundance (Y ^ 0.40), the envelope is consumed at a higher rate, because the 
hydrogen shell burns at higher temperatures (Dorman 1992; Sweigart & Gross 1976). The 
result is that more stars, with a greater range of envelope mass on the HB, will follow the 
PEAGB and AGBM paths (Fagotto et al. 1994c; DR093). It must be stressed that it is the 
high absolute value of the helium abundance that causes this effect (and not the increase in 
Z). 

High metallicity might affect the mass loss processes on the RGB. Gompared to 
metal-poor RGB stars, metal-rich RGB stars have lower masses and larger radii at the 
same luminosity, and therefore lower surface gravities, which may in turn increase the rate 
of mass loss on the RGB (Horch et al. 1992 and references therein). With greater mass loss, 
more stars will arrive on the EHB. We note that although the increase of mass loss with 
increasing metallicity is certainly plausible, the effect remains unsubstantiated by theory or 
observations. 

Once we understand the different physical parameters driving HB morphology, we may 
see how two different interpretations can arise for the correlation between Mg 2 strength 
and UV upturn strength. Arguing for the case of high metallicity, Bressan et al. (1994) 
say that more massive galaxies have higher mean and maximum metallicity, and thus — 
assuming Y also becomes high for high Z — more stars have a high enough metallicity to 
follow the evolutionary paths of the PEAGB and AGBM stars. Taking the opposing view, 
Lee (1994) argues that more massive galaxies have higher mean metallicity and are older 
than less massive galaxies, and thus whatever metal-poor fraction there is in the stellar 
population will have bluer HBs. Measurements of FUV line strengths have been proposed 
as sensitive tests of the models by both camps (see also Yi et al. 1995). In the Park & Lee 
(1997) model, one might expect to observe weak G IV and Si IV features anticorrelated 
with (mi 55 o — V). In the metal-rich scenario, one might expect strong lines and a positive 
correlation with (misso — V). 

The different metallicity scenarios have important implications for the ages of elliptical 
galaxies and in turn the age of the Universe. Bressan et al. (1994) demonstrated that in 
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the high-Z - high-F scheme, the onset of UV excess serves as an age probe, since it is the 
most rapidly evolving feature in the spectrum of an old elliptical galaxy. In their example, 
Bressan et ah (1994) show that the onset of UV excess occurs in the galaxy rest frame at 
an age of 7.6 x 10® yr. In the low-metallicity scheme. Park & Lee (1997) claim that giant 
elliptical galaxies must be at least 3 Gyr older than our own Galaxy, implying the lower 
limit to the age of the Universe must be 19 Gyr. Both models assume a Gaussian mass 
dispersion and either a Reimers (1975) mass loss law or hxed mass loss. The ages and 
metallicities at which the UV upturn begins to appear are unfortunately quite sensitive to 
these assumptions and to the parameters. The production of hot HB stars in the globular 
clusters and in the Galactic field is not easily described by a simple function, as we do not 
know how to characterize mass loss in stars. The UV upturn therefore has little value as an 
age indicator until the physics driving the phenomenon is understood. 

In contrast, Dorman et al. (1995) argue that the size of the population of EHB stars 
— i.e., the fraction of all HB stars in the galaxy that are UV-bright — is the only sure 
deduction that can be made from the FUV radiation. They estimate, by using broadband 
colors derived from simple models of all evolutionary phases, that the blue FUV sources 
must have EHB fractions reaching ~ 20 % of the HB population, a result only weakly 
dependent on the Z—abundance of the models. They argue that such a fraction is not 
a “trace population.” Instead, it is more plausible that the UV flux is produced by a 
signihcant minority of the dominant population, comprised of stars which have undergone 
higher mass loss than their red HB counterparts. The reason for the UV-Mg 2 correlation 
remains unknown, but may be explained by an increase in the likelihood of high mass loss 
in more metal-rich (or, at least, a— enhanced) stellar populations. In this picture there is 
no way to relate the FUV radiation to galaxy age, until mass loss in HB stars is better 
understood. In addition, the UV upturn is correlated with the velocity dispersion ((J„), and 
although the correlation between (misso — V) and Mg 2 index is tighter than that between 
(^^1550 — U) and (Burstein et al. 1988), it is not yet confirmed that the underlying 
correlation (i.e. causal relationship) is with abundance. 

Because of the great uncertainties in the physical mechanisms responsible for the 
production of EHB stars, we have chosen in this paper to look at the constraints provided 
by the UV spectra themselves, with only very broad constraints on the overall stellar 
populations in the galaxies. Our conclusions are thus independent of any particular model 
of elliptical galaxy evolution. The limitation of this approach is that without a model there 
is no straightforward way to specify the distribution of stars on the zero-age HB. This would 
be a serious problem if there were a wide variety of FUV spectral energy distributions 
for normal elliptical galaxies. However, observations indicate that E galaxy FUV spectra 
are remarkably uniform (Burstein et al. 1988; Brown, Ferguson, & Davidsen 1995) and 



consistent with a highly bimodal HB distribution (Ferguson 1994). Thus, for most purposes 
we can model the hot stellar population with only two evolutionary tracks, one for a PAGE 
star and one for an EHB star. We then adjust M^nv for the EHB star and the fraction of 
the population in each type of star to achieve the best £t. In §4.4 we discuss the effect of 
relaxing this assumption of extreme bimodality. 


1.2. Abundance Anomalies in Evolved Stars 

It is not a straightforward exercise to determine the abundances in the stellar 
populations producing the UV light from elliptical galaxies. This is because the UV 
absorption line strengths in the spectra of evolved stars do not provide a direct measure of 
their inherent abundances. Numerous observations of HB stars and their progeny in our 
own Galaxy show that these stars exhibit abundance anomalies (see Saffer & Liebert 1995 
and references therein). Helium is usually underabundant in sdB stars while overabundant 
in sdO stars. Garbon and silicon are usually underabundant in both sdB and sdO stars, 
while the nitrogen abundance usually appears normal in both classes (Saffer & Liebert 
1995; Lamontagne, Wesemael, & Fontaine 1987; Michaud, Vauclair, & Vauclair 1983). The 
FUV spectra of evolved stars observed with HUT show absorption lines that are markedly 
different from solar abundances or even scaled solar abundances (Brown, Ferguson, & 
Davidsen 1996). 

These abundance anomalies can be partly explained by models of diffusion processes in 
evolved stars, which demonstrate that abundances can be greatly enhanced or diminished 
through chemical separation in the stellar envelope and atmosphere. Individual elements 
can either levitate or sink in the outer layers of a star, depending upon which force has a 
stronger effect on each element: gravity (g) or radiative acceleration (grad)- In theory, the 
tendency is for heavy elements to sink in cooler stars, and for these elements to levitate 
in hotter stars. If grad is larger than g, an element will be radiatively supported and it 
accumulates in the atmosphere, increasing line saturation, which in turn will decrease grad, 
until an equilibrium can be reached between g and grad (Bergeron et ah 1988). If grad is less 
than g, the element will sink into the star. However, the resulting abundance decrease may 
or may not lead to an increase in the local value of grad, depending upon whether or not 
the lines were originally saturated or unsaturated. Thus, an equilibrium between g and grad 
may or may not be reached (Bergeron et ah 1988). Michaud et al. (1983) demonstrated 
that if the abundances of G, N, O, Ga, and Fe are all assumed to be originally at 1/200 
solar abundance, one can achieve abundance enhancements of three orders of magnitude 
through diffusion on the timescale of 10® yr. Unfortunately, diffusion theory has not reached 
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the level of accuracy needed to predict reliably the effects observed in evolved stars. For 
example, Bergeron et ah (1988) calculated three models with effective temperatures of 
20,000 K, 35,000 K, and 50,000 K, and with respective gravities of log g = 5.0, 5.5, and 
6.0, all originally at solar abundance. They found that while the nitrogen abundance was 
at approximately solar abundance for all models, agreeing with observations, the behavior 
of carbon and silicon did not agree with observed trends. In the models, carbon was 
underabundant and increasing with Tgfj, yet observations show that the carbon abundance 
seems to decrease with increasing Tefr (Bergeron et ah 1988 and references therein). Silicon 
was predicted to be between 0.1 and 1 of the solar value, but silicon dehciencies by factors 
of 10^-10^ are observed for these stars. Such discrepancies might point to other processes 
at work, such as a weak stellar wind (Bergeron et ah 1988). The bottom line is that these 
diffusion processes greatly complicate an absorption line analysis of an evolved population, 
and the strengths of the absorption lines are unlikely to reflect the inherent abundance of 
the population. Since the stronger metallic lines may well be saturated in stellar envelopes 
with very high abundances, metals in these envelopes may not be radiatively supported, 
allowing them to sink into the star. At the other extreme, very low abundances can be 
enhanced by orders of magnitude. Thus, diffusion processes tend to enhance low intrinsic 
abundances and diminish high intrinsic abundances. 

In the elliptical galaxies, line analysis is further complicated by the fact that the 
spectra represent composite systems with components that cover a range in temperature 
and gravity, which in turn affect line strengths just as abundance does. Furthermore, in 
the FUV, the density of lines is so large that only the very strongest lines will be clearly 
distinct in these composite spectra. 

With all of these complications in mind, an absorption line analysis alone will not 
answer the metallicity debate; we must look to the shape of the UV continuum as well. 
The shape of the continuum reflects the composite nature of these galaxies which no 
single temperature stellar model may duplicate (Brown et al. 1995), and thus reflects the 
evolutionary paths of the individual stars in the composite model. In this study, we use 
the stellar evolutionary models of DR093, Bressan et al. (1993), and Fagotto et al. (1994a, 
1994b, 1994c) to construct synthetic spectra for model stellar populations. The individual 
synthetic spectra which comprise the composite models are from the grids of Brown et 
al. (1996) and Kurucz (1992), for Tgff > 10,000 K and Tefr < 10,000 K respectively. The 
composite models which best £t the elliptical galaxy spectra draw the bulk of their flux 
from stars in the range 15,000 K < Teg < 30,000 K. 
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2. OBSERVATIONS 

In March of 1995, the Hopkins Ultraviolet Telescope (HUT) was flown aboard the 
Space Shuttle Endeavour during the Astro-2 mission. It collected FUV spectra of six 
elliptical and SO galaxies; these observations are summarized in Table 1. Together with 
the spectra of two galaxies from Astro-1, these data represent the only FUV spectra of 
early-type galaxies with wavelength coverage down to the Lyman limit at 912 A, and thus 
include the “turnover” in the spectra that is needed to fully characterize the spectral energy 
distribution. As determined by Burstein et al. (1988) with lUE, the (misso — V) color varies 
over a range 2.04-3.86 mag for these galaxies. All of them suffer from a minimal amount of 
reddening (0.00 < E{B — V) < 0.035 mag). 

HUT is uniquely suited for FUV spectroscopy of nearby galaxies. The telescope 
uses a 90 cm diameter //2 primary mirror and a prime focus, near-normal-incidence, 
Rowland-circle grating spectrograph to obtain spectrophotometry from 820 to 1840 A 
in first order with a resolution of 2-4 A. The detector consists of a microchannel plate 
coupled to a phosphor intensiher that is read by a Reticon into 2048 pixels. All six galaxies 
were observed during orbital night through a 10 x 56" slit. A complete description of the 
instrument and its calibration is given in Davidsen et al. (1992). Modifications and the 
resulting performance and calibration for the Astro-2 mission are described in Kruk et al. 
(1995). The final Astro-1 calibration of HUT is presented in Kruk et al. (1996). The fast 
focal ratio, large apertures (modihed somewhat for Astro-2), and wavelength coverage down 
to the Lyman limit make HUT an ideal instrument for observing extended objects and for 
determining the effective temperatures of hot UV sources. For Astro-2, an increase in the 
effective area by a factor of 2.3, combined with an extended duration flight, allowed HUT 
to collect much more data with higher signal-to-noise (S/N) than on Astro-1. 

Fig. 2 shows the HUT slit superimposed on two different images of the elliptical 
galaxy M 60 (NGC 4649) and the nearby spiral galaxy NGC 4647. The top panel is an 
optical image ( 5850-6820 A; Minkowski & Abell 1963) from the Digitized Sky Survey 
(DSS), and the bottom panel is an FUV image (1344-1698 A) from the Ultraviolet Imaging 
Telescope (courtesy of UIT; see also Ohl & O’Gonnell 1997). The figure demonstrates 
several important points. Of the six Astro-2 galaxies considered in this paper, M 60 is 
the only galaxy that has a possible close companion. However, the pairing is likely due to 
a superposition on the sky: although the two galaxies have similar redshifts, the lack of 
evidence of tidal interaction, either as tidal plumes or as a distortion of the spiral structure 
in NGG 4647, suggests that the galaxies are at different distances (Sandage & Bedke 1994). 
The bottom panel shows that the FUV light is very concentrated in M 60, and that a large 
percentage of this light falls within the HUT slit. The FUV image also shows that although 
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the FUV light from the spiral galaxy is more extended than that in the elliptical galaxy, 
our HUT spectrum is clearly not contaminated by the FUV flux from NGC 4647. 

Our six spectra were corrected for airglow contamination from the lines Lyman a A1216, 
Lyman (3 A1026, Lyman 7 A973, and O I A1305 by htting airglow prohles simultaneously 
with our model htting. The spectra were huxed and corrected for interstellar extinction 
in our Galaxy using the values of E{B — V) from Table 1, Ry = 3.05, and the extinction 
parameterization of Gardelli, Glayton, & Mathis (1989). For the stellar population 
model-htting described in this paper, the HUT data have been binned such that each 
wavelength point represents 5 HUT bins or 2.5 A. Only those portions of the observation 
with stable pointing and the galaxy well centered in the slit were used. Fig. 3 shows the 
HUT spectra for each of the six Astro-2 galaxies, along with models that will be described 
below. The NGG 1399 Astro-1 data has been shown in several previous papers (e.g. 
Ferguson et ah 1991). 


3. EVOLUTIONARY MODELS 

In order to model the FUV luminosity from the various classes of HB stars and 
their progeny, we construct synthetic composite spectra for evolved stellar populations 
by integrating individual synthetic stellar spectra over the stellar evolutionary tracks of 
DR093, Bressan et ah (1993), Fagotto et al. (1994a, 1994b, 1994c), and Vassiliadis & Wood 
(1994). The synthetic stellar spectra come from the grid of Brown et al. (1996) for effective 
temperature Teg > 10, 000 K, and from the grid of Kurucz (1992) for Teg < 10, 000 K 
(although these stars below 10,000 K do not signihcantly contribute to the integrated 
FUV luminosity). At each point along an evolutionary track, we hnd the synthetic stellar 
spectrum closest in Teg and surface gravity (g). To correct for the small mismatch in Teg 
and g, we scale the bolometric luminosity from the synthetic spectrum so that it matches 
that in the track step. We then weight the contribution from that point by the time 
spent at that location in the HR diagram. Our composite spectrum for a given track thus 
corresponds to a population that enters that track at the rate of one star per year. The 
normalization of the model and the distance to the galaxy then determine the actual rate 
of stars entering and leaving the evolutionary tracks, or the stellar evolutionary flux (SEF). 
Examples of these integrated composite spectra are shown in Fig. 1, for stars evolving along 
PAGB, PEAGB, and AGBM evolutionary paths. 

The diffusion processes in the outer layers of these evolved stars create inconsistencies 
between the abundances that determine the evolution of the stars and the abundances in the 
stellar atmospheres of these stars (see §1.2). Because of these inconsistencies, we consider 



12 


synthetic spectra with three different atmospheric metallicities {Z = Zq, Z = 0.1 Zq, 
and Z = 0.01 Zq) for all of the different evolutionary paths, regardless of the inherent 
abundances driving that evolutionary path. We define Zatm as the abundance of heavy 
elements driving the line strengths in our composite synthetic spectra, and we define Zevoi 
as the abundance of heavy elements driving the evolution of stars along the evolutionary 
path in a given composite synthetic spectrum. 

The limitations of our modeling scheme are the following: the small mismatch in 
temperature and surface gravity between the individual synthetic stellar spectra and the 
evolutionary track points; the assumption of local thermodynamic equilibrium in the 
synthetic spectra with Te^ < 45,000 K; pure H & He atmospheres for Teg > 20,000 K 
(with metals incorporated only in computing the emitted spectrum, not the atmospheric 
structures); the lack of continuum opacity for elements heavier than helium for spectra 
with Teff > 20,000 K; and the omission of mass loss during the HB and AGB phases in the 
DR093 evolutionary tracks, as well as the usual uncertainties and simplifying assumptions 
that go into stellar evolution calculations. Errors arising from the tracks are discussed in 
DR093 and in §4.3. These limitations may affect the quantitative results, but are less likely 
to affect the qualitative conclusions. 


4. MODEL FITTING 

As discussed in the previous section, the normalization of a model fit to the HUT data 
gives the stellar evolutionary flux (SEE) for stars evolving along the evolutionary path for 
that model (in the HUT slit). We also have a theoretical value of the SEFtot = the rate 
of stars evolving along all evolutionary paths in the galaxy. For old stellar populations 
(age > 10^° yr), the SEFtot per unit luminosity is approximately 2.2x10“^^ stars Lq 
yr“^ (Greggio & Renzini 1990). There is a dependence of this quantity on composition, 
particularly on the He abundance, since helium-rich stars at a given mass evolve more 
quickly. Our models assume that in a galaxy, the dominant population that controls the 
effective value of the SEFtot has Y close enough to Yq such that the solar value of the 
SEFtot is appropriate. To determine the bolometric luminosity for each of the Astro-2 
galaxies, we calculate the U-magnitudes through the HUT slit from archival HST WFPG 
images, correcting the 0.06 mag color offset between the F555W filter and the Johnson V 
band and using a bolometric correction (BG) of -1.32. This BG is consistent with a stellar 
population of age 12 Gyr and [Fe/H] = 0.25 (Worthey 1992). For a given evolutionary 
model fit to the HUT data, normalizing our measured SEF by the bolometric luminosity 
gives a distance-independent measurement of the fraction of stars evolving along that 
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evolutionary path, since we know the SEFtot- 

For each of the six Astro-2 galaxies, we perform two different classes of htting: 
FHB stars alone, or FHB and PAGB stars together. In all of our htting, we perform 
minimization with the IRAF routine SPFCFIT (Kriss 1994). For every model ht to the 
data, we assume that neutral hydrogen from the interstellar medium (ISM) absorbs at a 
column depth equivalent to that along a line of sight in onr Galaxy, as listed in Table 1. 


4.1. EHB tracks 

We initially ht FHB composite models based on evolntionary tracks from a variety 
of gronps (DR093; Bressan et ah 1993; Fagotto et ah 1994a, 1994b, 1994c). While 
FHB/AGBM seqnences from the latter two gronps were fonnd to prodnce reasonable 
hts to the HUT data, the DR093 tracks were fonnd to have somewhat snperior hts. In 
addition, they provide a hner mass spacing and the ability to check the ehects of Yzams and 
Zevoi separately (as well as the fact that we can compnte additional physically consistent 
seqnences where necessary). At hxed chemical composition and envelope mass, spectra 
constrncted with the Fagotto et ah (1994a, 1994b, 1994c) and Bressan et ah (1993) tracks 
agree with those constrncted from the DR093 tracks. Thns, to simplify the discussion, we 
limit onr analysis to the DR093 tracks. 

The DR093 tracks are broken into snbsets according to the abnndances of metals and 
helinm (see their Table 1 and onr Table 2); there are eight different abnndance combinations. 
The tracks detail the evolntion of HB stars from the zero-age horizontal branch throngh 
core He exhaustion, and continne either to the thermally pnlsing stage (for AGB evolntion) 
or to the white dwarf cooling track (for AGBM and PEAGB evolntion). Each track is 
parameterized by abnndance, ZAHB core mass, and ZAHB envelope mass. We took all 
136 tracks in their paper and constrncted a corresponding composite spectrnm for each 
track, which was then ht to the HUT data. These spectra were compnted with metallicities 
of Zatm = Zq, 0.1 Zq, 0.01 Zq. We fonnd that those tracks with Zatm = 0.1 Zq best 
matched the overall absorption featnres in the HUT data (see §6). The best ht evolntionary 
model (with Zatm = 0.1 Zq) for each of the Astro-2 galaxies is listed in Table 3 and plotted 
in Fig. 3. All of the best hts come from two track sets with high Z^voi and Yzams (sets G 
& H in Table 2). We also inclnde resnlts of the same htting applied to the spectrnm of 
NGG 1399, which was observed throngh a 9.4 x 116" slit on Astro-1 (we omit the M 31 
Astro-1 data because of the large and nncertain extinction toward that galaxy). The results 
in Table 3 show that the SEF for EHB stars (SEFehb) is a small fraction (^ 10 %) of the 
SEFtot- In a model comprised solely of EHB stars, one can reprodnce the amonnt of FUV 
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light from these galaxies without approaching the limits imposed by the fuel consumption 
theorem, since EHB stars and their progeny are such efficient FUV emitters. Using the 
statistics to judge the quality of ht, we can compare the quality of the best-htting EHB 
model from each abundance group (cf. Table 2) to that of the best ht of any group. We 
make this comparison in Table 3, which shows a strong preference for groups G and H in 
the htting. 

We shall use the FUV spectrum of M 60 to demonstrate how models with higher 
metahicity and helium abundance produce a closer match to the HUT data, since our M 60 
spectrum has the best S/N ratio. Fig. 4 shows the best ht Zatm = 0.1 Zq model for each 
of the 8 diherent abundance {Z^voi) sets. In the left column are the evolutionary models 
for stars with abundances less than or equal to the solar abundance (sets A-D in Table 2). 
These plots show a signihcant dehcit in the hux at the shortest wavelengths in the model, 
as compared to the HUT data for M 60. On the other hand, the evolutionary models in the 
right column, for stars evolving at higher metal and helium abundance (sets E-H in Table 
2), show a somewhat smaller flux dehcit at the shortest wavelengths. The evolutionary 
model that hts best is that with Z^voi = 0.04 and Yzams = 0.34 (set G). 

If we dehne the hux dehcit as 


AF 


970A 


i^data - F, 


model) 


X=912A 


970A 

^ ( Fdata 
X=912A 


we can characterize the dehcit numerically. These dehcits (listed in Fig. 4) clearly show that 
the low Yzams and Z^voi models in the left colnmn are more dehcient in this wavelength 
range, as compared to the high Yzams and Z^voi models in the right column. However, even 
the models in the right colnmn show a signihcant hux dehcit. This dehcit could be due 
to the absence of hux from PAGB stars in these models. Because of the high integrated 
Inminosity for stars evolving along EHB tracks (see Fig. 1), the SEFehb is a small fraction 
of the rate of all stars evolving along any evolutionary path in these galaxies (the SEFtot; 
cf. Table 3). Most of the stars must be evolving along the classical PAGB tracks. 

The low-metallicity models appear to have a stronger hux dehcit than the high- 
metahicity models; the reason may be related to the ratio of hux emitted in the HB and 
post-HB phases. We note that the mean temperatnre of the integrated spectrnm for a 
given evolutionary sequence depends on the ratio of the radiation emitted in the post-HB 
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stages to that emitted before core He exhaustion, simply because the mean temperature of 
AGBM stars is higher than the earlier stage. If one takes the integrated flux below 1200 A, 
and computes the ratio of this flux in the post-HB phase to that in the HB phase, one hnds 
that this ratio increases as Yzams and Z^voi increases. This is at least a partial explanation 
why our models tend to favor the higher metallicities. 

The models that work best from any of the abundance {Z^voi & Yzams) sets and for 
any Zatm are always those that integrate over AGBM evolutionary tracks (cf. Fig. 1). 
Models that integrate over PEAGB tracks have spectra that are too cool to agree with 
the HUT data. We hnd that a narrow range in envelope mass on the EHB produces 
spectra that agree reasonably well with the HUT data, consistent with the Endings in our 
previous analysis of these galaxies (Brown et ah 1995). This narrow distribution is on the 
blue end of the HB with low (but not the lowest) envelope masses on the HB. We discuss 
the distribution below in the context of the more realistic two-component EHB -|- PAGB 
models. 


4.2. EHB and PAGB tracks 

Because a very small fraction of the stars are evolving along EHB tracks in these 
galaxies, a more appropriate model to fit to the data is a two-component model, with 
both EHB stars and PAGB stars. Our PAGB models are constructed by integrating over 
the six solar-metallicity H-burning tracks of Vassiliadis & Wood (1994) with core masses 
=0.569, 0.597, 0.633, 0.677, 0.754, and 0.900 Mq. Vassiliadis & Wood (1994) also 
calculate H-burning tracks at lower metallicities (six at Z = 0.008, four at Z = 0.004, 
and two at Z = 0.001) and a handful of He-burning tracks. Since all of the H-burning 
tracks appear qualitatively similar, and since He-burning PAGB evolution is thought to 
be in the minority (^ 25%; Vassiliadis & Wood 1994 and references therein), we take the 
six H-burning solar-metallicity tracks as representative of the evolution for classical PAGB 
stars. We note that all of the Vassiliadis & Wood (1994) tracks start at Tefr = 10,000 K in 
the HR diagram, after the thermally pulsing (TP) stage (see Fig. 1). Our PAGB model 
spectra thus omit the flux that is generated as the star evolves along the AGB from the HB 
to the TP stage, since in the HUT wavelength range the contribution from the cool AGB is 
negligible. 

We used two methods to fit two-component (EHB & PAGB) models to the HUT data. 
The first method allowed the contributions of EHB stars and PAGB stars to float both 
freely in the fit. The second method constrained the contribution of the PAGB component 
so that it would not exceed the SEFtot- Assuming that the SEF^ot ~ 2.2 x 10“^^ stars 
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Lq yr“^ ±10% (Greggio & Renzini 1990), and that EHB stars by themselves can at most 
contribute a small fraction (^10%) of the SEFtot (from our own htting in the previous 
section, and from Brown et ah 1995), we restricted only the PAGB component of the £t, and 
not the sum of the EHB and PAGB components. Within the uncertainty of the SEFtot, the 
EHB contribution to the SEFtot is small. This assumption simplifies the minimization. 

When we did not restrict the PAGB contribution to the model fit, we found that EHB 
stars from all eight of the DR093 metallicity groups could reasonably fit the HUT data; 
the fitting only slightly favored the high Zevoi and high Yzams models (sets E-H in Table 
2). This was because the short-wavelength deficiencies of the models with Z^voi < Zq were 
compensated by a large PAGB contribution (^ half of the FUV flux). However, these 
solutions were not astrophysically consistent, since the PAGB contribution allowed in the 
fits was much larger than that allowed by the theoretical SEFtot (unless our models severely 
underestimate the potential UV flux from PAGB stars). 

When we restricted the PAGB contribution so that it could not exceed the theoretical 
SEFtot, we found that the EHB evolutionary models with high Z^voi and high Yzams were 
favored in the fits to the HUT data. The strength of this result is a function of the assumed 
in the galaxies, in the sense that the high Yzams and high Zevoi models are more 
strongly favored as the assumed increases. Our results depend upon 

because as the mass increases, the PAGB stars evolve more rapidly and produce less FUV 
flux over their lifetimes. Thus, if the stellar evolutionary flux of PAGB stars (SEFpagb) for 
each is limited by the same SEFtot, the more massive PAGB stars will contribute 

less flux to the EHB ± PAGB models than less massive PAGB stars. 

Since the strength of our results depends upon the assumed, we show these 

results for the two lowest masses in the set of PAGB models: those with = 0.569 Mq 

(Table 4 and Fig. 5) and 0.597 Mq (Table 5). Our results are stronger as increases, 

and thus the results in Table 4 represent the conservative results, while those in Table 5 
use a that is closer to that modelled in the planetary nebula luminosity functions 

of ellipticals (0.61 Mq; see below). In these EHB ± PAGB models, the SEFehb is a small 
fraction of the SEFtot predicted by the fuel consumption theorem, while the SEFpagb is 
much larger. Note that consistent models only require that SEFpagb < SEFtot, and not 
SEFehb ± SEFpagb = SEFtot- The value of SEFpagb that we derive depends on the 
integrated flux of a particular assumed PAGB track. In reality the PAGB population will 
have a mass spectrum whose mean mass can be higher and it integrated flux therefore 
lower, in which case our low-mass PAGB models will underestimate the number of PAGB 
stars present while overestimating the mean UV contribution from individual stars. The 
individual PAGB tracks have very similar integrated spectra, and so we cannot constrain 
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the PAGE mass with our fitting. 

Using the statistics to judge the quality of ht, we compare in Tables 4 and 5 the 
quality of the best-htting EHB + PAGE models from each abundance group (cf. Table 2) 
to that of the best £t. In Table 4, evolutionary abundances of at least solar are favored 
in 5 of the 7 galaxies. The exceptions are NGG 3115 and NGG 3379, where there are no 
preferences for either high or low evolutionary abundances in the fitting. Alternatively, in 
Table 5, evolutionary abundances greater than the solar value are favored in all 7 of the 
galaxies, with most fits favoring high Zf,„oi and Yevoi- When was increased to even 

higher values, the fits for all of the galaxies increasingly favored higher Ze„oz and Yzams- 
The culmination of this trend was shown in the previous section, where we fit models 
consisting solely of EHE stars. Such models are equivalent to the assumption that the 
PAGE stars are too massive to contribute significant flux to the FUV spectra. 

What is the appropriate to use in our fitting? Studies of the planetary 

nebula luminosity function (PNLF) may indicate that ^ 0.61 Mq. The PNLF is 

well-studied in both spiral bulges and elliptical galaxies. Giardullo et al. (1989) calibrated 
the PNLF in the bulge of M 31, and found that it was in excellent agreement with a 
set of PAGE tracks having a mean of 0.61 Mq. Since then, the same PNLF 

has demonstrated very good agreement with that in the bulge of the Sb galaxy M 81 
(Jacoby et al. 1989), in E and SO galaxies of the Leo I Group (Giardullo, Jacoby, & Ford 
et al. 1989), and in E and SO galaxies of the Virgo Gluster (Jacoby, Giardullo, & Ford 
1990). Four of the Astro-2 galaxies (M 49, M 60, M 87, and NGG 3379) were included in 
these PNLF studies. In the inner regions (isophotal radius < 1') of the Virgo and Leo I 
ellipticals, PN were lost in the bright galaxy background, but no significant gradients in the 
luminosity-specific PN density were seen beyond these regions. Thus, we can not know for 
certain the characteristics of the PAGE population within the cores of the ellipticals, but 
the evidence we have points to 0.61 Mq. 

In the literature there are other sets of PAGE tracks that extend the range of core 
masses even lower than the 0.569 Mq track of Vassiliadis & Wood (1994) and thus increase 
the potential PAGE contribution to the total UV flux. Although it is unlikely that all of 
the PAGE stars in a galaxy would be evolving along such low-mass tracks (see above), we 
investigated the effect this would have on our modeling. For example, the Schonberner 
(1987) set has a track with a core mass of 0.546 Mq. Such a low-mass PAGE star evolves 
so slowly that it contributes significant flux as the PAGE component of an EHE -|- PAGE 
model, even when the Schbberner track is restricted by the SEFtot (although we stress that 
the integrated flux from this 0.546 Mq track does not £t the HUT data by itself, without 
an EHB/AGBM component). So, EHB -|- PAGE models that used a 0.546 Mq PAGE 
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track did not show a preference for either the high Z^voi or low Z^voi models; the results were 
similar to those of the unconstrained EHB + PAGE models above. However, Schonberner 
states that this 0.546 Mq model is below the threshold for the occurrence of thermal pulses; 
the model is therefore actually a PEACE sequence. Although pre-TP evolution contributes 
negligibly to the integrated FUV luminosity of stars evolving along PAGE tracks, the same 
cannot be said of all PEACE stars. If we include the luminosity from a 0.550 Mq AGE 
track (DR093) as a plausible example of the pre-TP evolution for a star evolving along 
the 0.546 Mq Schonberner track, we hud that the composite spectrum appears similar to 
the PAGE spectrum in Fig. 1, except that the spectrum rises considerably at wavelengths 
longer than 1500 A, and thus does not agree with the Astro-2 data. The result is that 
this “PAGE” contribution (which is really PEACE) in an EHE -|- PAGE £t to the data 
is minimized. In summary, the 0.546 Mq track is an acceptable PAGE component only 
if two conditions are met. First, all of the PAGE stars in the galaxies must evolve along 
extremely low-mass tracks. Second, the pre-TP evolution must be characterized by high 
mass stars that are too cool to produce flux in the HUT range while on the AGE, yet lose 
enough mass on the AGE to subsequently become 0.546 Mq PAGE stars. 

Because the “pseudo-continuum” shape is affected by Zatmi the best htting model was 
usually, but not always, that with Zatm = 0.1 Zq. However, at any Zatm) the fitting still 
favored those evolutionary models with high Zevoi and high Yzams- Since the individual 
absorption features are better matched by Zatm = 0.1 Zq (see §6), we chose to demonstrate 
the best fits to the HUT data by consistently using models with this atmospheric abundance. 

Although the models in Fig. 5 fit the Astro-2 data fairly well, there are a few curious 
discrepancies. The M 49 data shows a large excess at 961 A, although the deviation from 
the model is not very significant. The formal probability of such a deviation from the 
model flux is 0.16. However, a more significant discrepancy appears in the NGG 3379 
data as an emission feature near G IVAA1548,1551. The feature can be fit by a Gaussian 
with A = 1562.7 A, equivalent width = 10 A, and FWHM = 3000 km The statistical 
significance of the feature is 0.0084 (better than a 2a detection). But what could be 
causing this feature? We observed NGG 3379 twice, and the same feature shows up in both 
observations, so it is unlikely that some “glitch” in the data is causing the feature. The 
deredshifted wavelength of the feature would be 1557.9 A, which is still 7 A longward of 
the most likely emission at G IV. If an extraneous emission source fell within the HUT slit, 
such as a supernova remnant or cooling flow, one would expect to find other emission lines 
(Lyo;, O VI, etc.) at the same redshift. Given the lack of other detectable emission features 
in the NGG 3379 spectrum, we tentatively conclude the feature is a fluke; a 2a fluke is not 
unreasonable given 6 data sets having 410 bins each. 
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The best-fitting two-component models incorporate EHB tracks with a narrow range in 
temperature and envelope mass on the EHB, but this range does not include the very lowest 
envelope masses, which extend down to 0.002 Mq. Since the contours are probably 
dominated by the systematic uncertainties of the models, we conservatively use the 
statistics to determine the range of M^nv that is within 4a of the best-£t model. We find 
that in our two-component EHB -|- PAGB models, a narrow distribution of M^nv on the 
EHB fits the Astro-2 data within 4a, as we found in our initial analysis of these data using 
lower resolution (10 A) models (Brown et al. 1995). The range for each galaxy is listed in 
Tables 4 and 5. 

As in the previous section with the one-component (EHB) models, we use the high 
S/N data for M 60 to demonstrate why models with high Zevoi and high Yzams work better 
than models with low Zevoi and low Yzams- Fig- 6 shows the best fit EHB -|- PAGB model 
from Table 4, for each of the 8 different abundance {Zevoi) sets. In each panel, the PAGB 
component has Mcore = 0.569 Mq. The flux deficiency (AE) at the shortest wavelengths 
is more significant in the low Zevoi and low Yzams models than in the high Zevoi and high 
Yzams niodels. However, the high Zevoi & Yzams niodels do still show some deficiency at 
the shortest wavelengths, and we discuss this problem in the next section. 

In our previous study of these galaxies (Brown et al. 1995), we found that PAGB stars 
alone do not fit the FUV spectral energy distribution of any of the galaxies. This result also 
holds true in the current study. The result that none of the galaxies can be fit by PAGB 
stars is significant, because if one were trying to reproduce only the flux at 1550 A, and 
not the entire FUV spectral energy distribution observed with HUT, galaxies with very red 
(^^1550 ~ y) colors could be accounted for by PAGB stars alone. Hence all of the galaxies 
in our study appear to require some EHB fraction. 


4.3. Systematic Errors 

As we discussed above, HB models from DR093 with supersolar Yzams and Zevoi ht 
the HUT data better than models with subsolar Yzams and Zevoi- This is true for both 
EHB models and EHB -|- PAGB models. Although the supersolar Yzams and Zevoi models 
work better across the entire HUT range, the largest improvement is in the region shortward 
of 970 A. But how much can systematic errors contribute to this result? The major sources 
of systematic error in our modeling are errors in the individual stellar synthetic spectra, 
uncertainty in the extinction correction for the galaxies, and uncertainties in the behavior 
of the input stellar evolution models. 
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As described in Brown et al. (1996), the Lyman series lines are somewhat too strong 
in our grid of synthetic stellar spectra. This produces a dehciency in the synthetic spectra 
as one moves to shorter wavelengths in the HUT range. In Brown et al. (1996), there are 
six plotted synthetic stellar spectra that have flux down to the Lyman limit, at effective 
temperatures of 17000 K, 24000 K, 29900 K, 36100 K, 40000 K, and 55000 K, and each has 
HUT data or a Kurucz LTE model for comparison. If we measure the flux dehciency AF for 
these 6 models, we hnd that the average dehciency is 0.12±0.08. Therefore, the supersolar 
Yzams & Zevoi models in the right column of Fig. 6 have dehciencies that fall within the 
range expected from systematic errors in our individual synthetic spectra, whereas the 
subsolar Yzams & Zevoi models in the left column have dehciencies that fall outside of that 
range. 

The extinction parameterization of Cardelli et al. (1989) does not extend below 1000 A, 
but we have extrapolated it to the Lyman limit at 912 A. Although this results in some 
uncertainty in the extinction correction, we note that our extinctions are quite low (see 
Table 1), and thus any errors in the shape of the extinction curve will in turn be quite 
small. For example, our use of Ry = 3.05, characteristic of the very dihuse ISM, produces 
an extinction correction that deviates from the galactic mean curve {Ry = 3.1) by less than 
1% when E{B — V) = 0.035, which is the highest extinction in Table 1. The significant 
uncertainty in the extinction is the possibility of extinction intrinsic to the galaxies. 
However, we note that if the Astro-2 galaxies have internal extinction, our correction for 
this extinction would increase the flux deficits for all of the models, and thus our supersolar 
Yzams & Zevoi models would still match the Astro-2 data better than those models with 
subsolar Yzams & Zevoi- 

The DRO90 tracks neglect the effects of mass loss. If mass loss takes place by winds 
during the EHB evolution, it is likely to affect the brighter, low-gravity AGBM stage. It 
will drive the evolution to higher temperatures faster, but since it reduces the available 
hydrogen fuel it is also likely to reduce the total UV output in these stages. The neglect of 
mass loss in the tracks should not affect our conclusions regarding metallicity. 


4.4. Mass Distribution on the HB 

The composite spectrum for a given evolutionary track is the flux one would observe 
from a population of stars evolving from a single mass on the HB. In reality, the FUV 
spectra in elliptical galaxies probably arise from populations on the HB (and their progeny) 
with some distribution of mass. Unfortunately, the true mass distribution on the HB 
and the mechanisms that govern it remain highly uncertain. These uncertainties are 
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compounded by the uncertainties in the age and composition of these stellar populations. 
Thus, there is scant theoretical justihcation for any given mass distribution. 

As an initial attempt to constrain this problem, we tried to £t Gaussian distributions 
of envelope mass (using a single core mass) to construct our populations, but the best 
fit always occurred for a ^-function in envelope mass, with the quality of the fits rapidly 
declining as the width of the Gaussian increased. We believed that a grid of evolutionary 
tracks with a finer distribution in envelope mass would allow a Gaussian distribution to fit, 
so we reran the evolutionary code of DR093 to create a new subset of tracks. These 11 
AGBM tracks were evenly spaced by AM^nv = 0.002 Mq around the best fit EHB track for 
M 60 in Table 3 {M^nv = 0.036 Mq). 

We then performed another two-component EHB -|- PAGB fit to the M 60 HUT data, 
using the same Zatm and Z^voi as the best EHB -|- PAGB fit in Table 5, but now the EHB 
contribution came from a Gaussian distribution of envelope mass on the HB, such that 
the number of stars entering bins of unit envelope mass were constrained by a symmetric 
Gaussian. We found that a narrow distribution of M^nv, centered on M^nv = 0.042 Mq 
with FWHM = 0.003 Mq, fit the M 60 data slightly better than the original ^-function 
at Menv = 0.046 Mq. While this result suggests that the HB distribution is extremely 
bimodal, the systematic uncertainties in the model atmospheres and the evolutionary tracks 
combine to make the constraints weaker than they appear from formal fitting. Model 
spectra constrained with a fairly flat HB distribution are not very different from those with 
a pure AGBM component, because the long lifetimes make the AGBM stars dominant. 
For example, increasing the envelope mass FWHM from 0.003 Mq to 0.1 Mq increases 
by 20, and increasing the FWHM to 0.5 Mq increases x^ by 32, but the SED in these 
models still provides a reasonable match to the continuum shape of the HUT data. Thus 
a “flatter” distribution in mass can produce a reasonable SED because the distribution in 
flux contribution is still highly peaked around the AGBM stars. The Aa limits on the Menv 
distribution given in Tables 5 and 6 denote the limits on which stars can dominate the SED. 


5. COLOR-COLOR ANALYSIS 

The (mi55o — V) color of elliptical galaxies is the standard measure of the strength of 
the UV upturn. If galaxies with stronger UV upturns have a larger fraction of their stellar 
populations evolving along EHB evolutionary paths than galaxies with weaker UV upturns, 
it then follows that galaxies with stronger UV upturns should also have a smaller fraction 
of stars evolving along the classical PAGB tracks (see Brown et al. 1995 and Ferguson & 
Davidsen 1993). While the (misso — V) color tracks the strength of the FUV flux relative to 
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the optical flux, the (m<iooo ~ ''^1550) color tracks the “characteristic” temperature of the 
FUV population. Using the HUT data and archival HST images, we define these colors as 

(^<1000 ~ ''^1550) = 2.5 logic — 

•^A=912-1000A 

and 

(mi55o -V) = 2.5 logic —o - — -7^ 

■^A=1450-1650A 

where log fvo=io = -12.41 ergs cm“^ s“^ (Johnson 1966), and Vo is derived from the V 
and E{B — V) values in Table 1. 

The left-hand panel of Fig. 7 shows the color-color diagram for the six Astro-2 galaxies 
and for NGC 1399 (we again omit the Astro-1 galaxy M 31 because of the large and 
uncertain extinction toward that galaxy). The dashed line shows the shift in (misso — V) 
that would occur if we used the measurements of Burstein al. (1988) through the smaller 
lUE apertures. The large extent and variation in the dashed lines is a clear indication of 
the strong variation in the FUV/optical flux ratio as a function of radius in these galaxies. 
These differences in color gradients are apparent in the simultaneous UIT observations of 
these galaxies (O’Connell et al. 1992; Ohl & O’Connell 1997). 

In the right-hand panel, we show the variation in (m<iooo — for the integrated 

spectra of models in abundance group C (Table 2). The four reddest of these models 
(DR093) stop at the TP stage (and produce negligible FUV flux), so we have tacked on 
plausible PACE evolution from Vassiliadis & Wood (1994). The (m<iooo — ^ 1550 ) color 
appears “hot” for the high and low extremes of HB mass (respectively HB/PACB and 
EHB/ACBM). As one moves from the blue end of the HB to larger envelope masses, the 
flux from the EHB stars and their progeny becomes progressively “cooler” in the FUV. 
This trend may explain the positions of the HUT data points in the left-hand panel. The 
hot NCC 1399 spectrum may be dominated by stars evolving from the blue extreme of the 
HB, while the hot FUV spectra of M 49, NCC 3115, and NCC 3379 may have significant 
contribution from PACE stars. The “cooler” FUV spectra of M 60 and M 89 may be 
dominated by ACBM stars with higher envelope masses than those in NCC 1399. Of 
course, with only a handful of galaxies, it is premature to deduce a correlation between the 
behavior seen in the data and that shown in the right-hand panel of Fig. 7. But this type 
of color-color diagram might be a useful tool in future studies. 

M 87, which has an active nucleus, might have other reasons for its position in the 
color-color diagram. The active galaxies in the Burstein et al. (1988) study did not follow 
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the UV-Mg 2 relation of the quiescent galaxies. However, recent HST observations have 
shown that the FUV flux from the point source and jet in the M 87 nucleus cannot be 
contributing signihcantly to the HUT spectrum (Tsvetanov et ah 1996). Also, subtracting 
the nuclear flux from the V measurement through the HUT slit would only diminish V and 
make the (misso — V) color even bluer. So why is M 87 so blue in both of the colors in 
Fig. 7? We believe that the answer may be related to the position of M 87 in the Virgo 
cluster. Because M 87 lies at the center of this cluster, it has been accreting gas and 
metals from other galaxies. Perhaps this accretion has altered the stellar evolution in M 87 
so that it is markedly different from that in the other Astro-2 galaxies. For example, a 
recent episode of star formation could dramatically decrease (m<iooo — but not affect 

(^^1550 ~ V) signihcantly. If we start with the best ht EHB + PAGB model for M 60, and 
add a third component representing continuous star formation (SF) for the past 10^ yr at 
the rate of 0.16 Mq yr“^ with a normal (Salpeter) IMF, we hnd that we can reproduce the 
M 87 HUT hux. Although the total accretion rate from the M 87 cooling how is estimated 
at 20-30 Mq yr“^ (White & Sarazin 1988), the rate of star formation is generally a small 
fraction of the cooling rate in giant ellipticals; the fraction has been calculated as high as 
15 % in some cases, but 1 % in many others (O’Connell & McNamara 1988). In our 
EHB + PAGB + SF model, SEFehb = 2.8x10“^ stars yr“^ and SEFpagb = 0.36 stars 
yr“^. If we compare these numbers to those in Tables 4 and 5, we see that adding the 
young stars means the ratio of EHB to PAGB stars in M 87 would be more in line with the 
ratio in M 60. Although this small number of young stars would make (m<iooo ~ ''^ 1550 ) 
signihcantly bluer, their contribution to the V band would be negligible, since this hux is 
dominated by a much larger number of RGB stars. Unfortunately, the S/N in the M 87 
data is low enough that we cannot put limits on the contribution of young stars using the 
C IV and Si IV absorption. In contrast, the NGC 1399 data is of higher S/N, and the lack 
of detectable C IV and Si IV absorption in the NGC 1399 spectrum makes star formation 
an unlikely explanation for its FUV colors (cf. Ferguson et al. 1991). 


6. LINE STRENGTHS 

Although the line strengths in our FUV spectra probably do not rehect the inherent 
abundances driving the evolution in these galaxies, they can still tell us information about 
the abundances in the atmospheres of stars in these evolved populations. We found that 
EHB + PAGB models with Zatm = 0.1 Zq best reproduced the overall absorption features 
in the HUT data. Fig. 8 shows the M 60 data and the best fit EHB + PAGB model for 
M 60 (Table 4), but with three different atmospheric metallicities: Zatm = Zq, 0.1 Zq, 
and 0.01 Zq. It is obvious that the practically featureless 0.01 Zq model (bottom panel) 
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has absorption lines that are too weak to match the HUT data. In contrast, the Zq model 
(top panel) has absorption lines which are mostly too strong. Most of the absorption lines 
are best matched in the 0.1 Zq model (center panel), although the broad blanketing from 
Fe absorption at 1550 A and 1000 A might be better matched in the Zq model. Thus, in 
M 60, the atmospheric metal abundance is likely to be somewhere in the range 0.1-1 Zq. 

To quantify the atmospheric abundances in these populations and to look for any 
correlation between the line strengths and colors, we measure the absorption line strengths 
empirically, by absorption line htting and spectral indices. These two methods do not 
rely upon stellar evolutionary models or synthetic spectra. However, in practice, the low 
S/N in the HUT data hampers these measurements and makes it difficult to draw strong 
conclusions regarding the atmospheric abundances in these populations. 


6.1. Spectral Indices 

By using the EHB models as a guide, we can dehne spectral indices which reflect the 
abundances of elements responsible for strong absorption features in the FUV: C, N, Si, and 
Fe. While C, N, and Si each have a small number of strong features in the HUT wavelength 
range, Fe produces strong line blanketing over large regions of this range (Brown et ah 
1996). 

We dehne line and continuum regions for C, N, Si, and Fe, and then measure the 
spectral index I in magnitudes as 
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where Ncont and Nums are the number of continuum bins and line bins, respectively. 

Table 6 lists the line and continuum regions used to dehne our spectral indices. Because 
the FUV is completely blanketed by absorption lines, there is no way to dehne a “clean” 
spectral index that will be ahected by absorption from only a single element (see Figs. 9-12 
in Brown et ah 1996). To show how these indices may change as a function of abundance 
in a stellar atmosphere, we measured the indices in an LTE synthetic stellar spectrum 
with Teff = 25,000 K, log g = 5.0, and Zatm = Zq. We repeated this measurement in 
synthetic spectra where one of four elements (C, N, Si, and Fe) was missing, and the results 
are plotted in Fig. 9. Because the indices are not completely “clean,” the index from one 
element can be ahected by the abundance of other elements (see for example the ehect of 
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removing Fe in the upper right panel of Fig. 9). However, it is apparent in the hgure that 
these indices are sensitive to the individual abundance changes in these four elements. 

To improve the S/N, we combined the deredshifted spectra from galaxies with strong 
UV upturns (M 60 and M 89) to produce a “UV-strong” spectrum, and the deredshifted 
spectra from galaxies with weak UV upturns (M 49, NGC 3115, and NGC 3379) to produce 
a “UV-weak” spectrum. The individual spectra were weighted by the variance in this 
process. We then measured the spectral indices in both the UV-strong and UV-weak 
spectra. The UV-strong and UV-weak spectra are plotted in Fig. 10, along with the regions 
used for the spectral index measurements. 

Fig. 11 shows the variation in these indices vs. Zatmi for each of the EHB -|- PAGE 
models shown in Fig. 6 and Table 4. In this way, we can account for the spread in indices 
due to the uncertainties in Zevoi- each panel shows the indices as a function of Zatm for a 
model from each of the Z^voi groups in Table 2. On the left-hand side of each panel, we plot 
the indices for the UV-weak and UV-strong Astro-2 data. These indices are shown as error 
bars that encompass the la statistical uncertainties. Because these uncertainties are fairly 
large, we cannot use the indices to determine accurately what the atmospheric abundances 
are in the stellar populations in these galaxies. However, it is evident in Fig. 11 that the 
indices are consistent with Zatm ~ 0.1 Zq. It is also apparent that the HUT data exhibit 
a weak trend, such that the UV-strong galaxies tend to have stronger UV spectral indices 
than the UV-weak galaxies; seven of the nine indices in Fig. 11 show this tendency. 


6.2. Absorption Line Fitting 

We attempted to £t absorption lines in the two Astro-2 galaxies with reasonable S/N: 
M 60 and NGG 3379. Each of these galaxies was observed multiple times during the 
Astro-2 mission, in order to improve the S/N in two galaxies with very different UV upturn 
strengths. We also fit these same features in the best-fit EHB -|- PAGE model for M 60 
{Zatm = 0.1 Zq] Table 4). The fits to the two galaxies and the model are listed in Table 7. 

The low S/N of the Astro-2 data required that we restrict the number of free parameters 
in our absorption line fits. Thus, we fit Gaussian absorption profiles to those features listed 
in Table 7, holding the FWHM and the line center fixed. The FWHM was held to the 
quadrature sum of the intrinsic line width (determined from the EHB -|- PAGE model), 
the velocity dispersion in the nucleus of each galaxy, and the 3 A resolution of the HUT 
data. The velocity dispersions of M 60 and NGG 3379 are 363 and 214 km s“^ respectively 
(Burstein et ah 1988). The uncertainties in equivalent width (EW) are calculated from 
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the la errors, the upper limits in EW are calculated for the 2a conhdence level, and the 
signihcance is the formal probability that the absorption feature results from a statistical 
fluctuation. 

An analysis of the line strengths in Table 7 does not provide any indication of 
systematic changes between the absorption lines in the UV-strong M 60 and those in the 
UV-weak NGC 3379. However, given the statistical uncertainty in these measurements, a 
correlation between the FUV colors and the FUV absorption line strengths could easily be 
lost in the noise. The absorption lines are again consistent with Zatm ~ 0.1 Zq. 


7. DISCUSSION AND SUMMARY 

Each individual analytical technique we have employed provides evidence for somewhat 
weak conclusions about the chemical composition of the stellar populations in these galaxies. 
Taken in their entirety, our results favor the hypothesis that the FUV light from elliptical 
galaxies originates in a population with high metallicity and possibly high He abundance. 
These conclusions appear to depend on the ratio of post-HB to EHB flux derived from 
the tracks. The results for the tracks with Z^voi > Zq are less distinguished from each 
other than they are from those of lower metallicity in this respect, so that more specific 
conclusions cannot be drawn. We cannot, however, from the FUV spectra alone, rule out 
the low-metallicity scenario of Park & Lee (1997). 

The preceding sections have demonstrated: 

1) The HUT Astro-1 and Astro-2 E galaxy spectra are best ht by a composite EHB -|- PAGB 
star model; the FUV spectra are dominated by EHB stars following AGBM evolution. 

2) When the flux from the PAGB component is constrained by the fuel consumption 
theorem, EHB stars of high Yzams and Z^voi are favored to make up the flux dehcit near 
the Lyman limit. 

3) Absorption features are consistent with Zatm — 0.1 Zq and may tend to increase as the 
UV upturn increases, but the direct constraints on metallicity from spectral lines are weak 
due to the low S/N of the data and the likely redistribution of elements within the EHB 
star atmospheres. 

If EHB stars alone could account for the FUV light in elliptical galaxies, then EHB 
stars with high Zevoi and high Yzams would best fit the Astro-2 data. However, we know 
that PAGB stars exist in ellipticals; we can observe bright PAGB stars in elliptical galaxies, 
and the fuel consumption theorem tells us that the stars evolving from the EHB are only 
a small fraction of the entire evolved population. If all of the stars in the population 
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evolved from the EHB, the FUV light from ellipticals would be much brighter than that 
observed, since EHB stars are such efficient FUV emitters. Dorman et ah (1995) estimate 
(^1550 ~ y) 0 for systems where all of the HB stars lose sufficient mass on the RGB 
to become EHB and subsequently AGBM stars (i.e. sdB and then sdO stars). Hence the 
theoretical bluest possible (misso — V) color implies ~ 6 times more FUV radiation than is 
observed in the nncleus of blnest galaxy, NGG 1399. Since we know PAGB stars make np 
the bnlk of the FUV population and contribute to at least a fraction of the FUV light, a 
more realistic approach is a two-component EHB -|- PAGB model. If the EHB and PAGB 
components are fit to the HUT data without fuel consumption constraints, then the HUT 
data can be reproduced equally well by populations with low Zg^oi & Yzams or high Z^voi & 
Yzams- The ambiguity occurs because a “hot” PAGB spectrum can offset the ffux dehcit 
at short wavelengths in the low Z^^oi and Yzams models. However, an nnconstrained EHB 
-|- PAGB model is unrealistic, because in snch hts the number of PAGB stars exceeds that 
allowed by the fuel consumption theorem. Although the EHB and nnconstrained EHB -|- 
PAGB models are instrnctive, a more realistic model is one where the nnmber of PAGB 
stars is astrophysically plausible and thus constrained by the fuel consumption theorem. 
When we £t these constrained EHB -|- PAGB models to the HUT data, the fitting favors a 
population evolving with high Zevoi and high Yzams- 

The EHB stars that best fit the Astro-1 and Astro-2 data in any EHB or EHB -|- PAGB 
model are those that lie within a narrow distribntion of envelope mass on the low-mass 
(blue) end of the HB. The stars in this narrow distribntion evolve along AGBM paths, 
bnt the mass distribntion does not include those EHB stars with the very lowest envelope 
masses on the HB. In addition, since the SEDs of PEAGB stars are so dissimilar to the 
Astro-2 data (see Fig. 1), PEAGB stars cannot comprise a significant fraction of the stellar 
popnlations in these galaxies. This result is partly to be expected from the small mass 
range that produces them for a given Y and Z. However, it is possible (and plausible) that 
relatively low mass PEAGB stars (with log L/Lq ~ 3.2) are produced in old populations 
through mass loss on the AGB. Such objects are in principle capable of snpplying the total 
UV flux for the weaker UV upturn systems (with misoo — V 3.6). Our modeling implies 
that all of the galaxies we have stndied instead have some significant contribntion from a 
small EHB population. From this we infer that the galaxies have a bimodal HB population 
in general. Nesci & Perola (1985) showed that a strongly bimodal temperatnre distribution 
on the HB was needed to explain the 2200 A dip seen in lUE spectra of elliptical galaxies, 
and Ferguson (1994) invoked a bimodal temperatnre distribntion to explain the SEDs of 
the Astro-1 observations of NGG 1399 and M 31, since a uniform distribntion of mass on 
the HB would produce a “flatter” spectrum (Ferguson 1994) than observed. 

The open cluster NGG 6791 might be an analog of the sitnation in elliptical galaxies. 



since it shows a high-metallicity evolved population with a strongly bimodal temperature 
and mass distribution on the HB (Liebert et al. 1994). However, the existence of this cluster 
might present problems to the theorists in both the high-metallicity and low-metallicity 
camps in the elliptical galaxy composition debate. It is also possible that this cluster 
may not be providing us direct information about mass loss in the single stars of metal 
rich systems, since several of the observed stars are binaries (E. M. Green 1996, private 
communication) which may have affected the evolution of the giant before core helium 
ignition, by tidal interaction and by mass ejection. However, bimodality occurs in globular 
cluster HB morphology as well (e.g. NGC 1851; Walker 1992), where binarism is much 
less common. We currently lack the understanding to explain the existence of a bimodal 
mass distribution in a coeval, single-metallicity population. The mass-loss mechanisms 
on the RGB that determine the distribution of mass on the HB are the largest source of 
uncertainty in the evolution of these populations. 

Given a bimodal HB distribution where the efficient UV emitters are on the EHB, 
the (mi 55 o — V) color should track that fraction of stars evolving from the EHB, while the 
(m<iooo — ''^ 1550 ) color should track the “characteristic” temperature of the FUV population. 
The uncertainties in mass distribution on the HB means that we will need more than 
seven data points to understand the behavior in a color-color diagram like that shown in 
Fig. 7. However, these FUV colors are consistent with the scenario proposed by Ferguson 
& Davidsen (1993) that UV-weak galaxies have a relatively large PAGB contribution to 
their spectra, while the UV-strong galaxies are dominated by EHB stars and their AGBM 
progeny, evolving from a range of temperatures on the blue end of the HB. In the future, 
with more data points in such a color-color diagram, we might be able to deduce information 
about the HB mass distribution. 

Park & Lee (1997) argue that stars in the low-metallicity tail of an extended metallicity 
distribution are responsible for the UV upturn seen in elliptical galaxies. Bressan et al. 
(1994) believe that a high-metallicity population produces the UV flux in these galaxies, 
and they explore metallicities as high as Z = 0.05. If the absorption features in evolved 
populations reflect their intrinsic abundances, then the actual situation in elliptical galaxies 
lies somewhere in between these two extremes, since the absorption features in the HUT 
data are consistent with Zatm — 0.1 Zq. There is a weak trend in the spectra of these 
galaxies, such that galaxies with stronger UV upturns show stronger UV line indices. This 
trend is in the opposite sense to that predicted by Park & Lee (1997), but in agreement 
with the predictions of Bressan et al. (1994) and Dorman et al. (1995). 

However, diffusion processes in the envelopes of stars in evolved populations tend to 
enhance low intrinsic abundances and diminish high intrinsic abundances, and thus it is 
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unlikely that Zatm = Z^voi- Given this situation, we must also compare the spectral energy 
distributions predicted by modeling the FUV data. This comparison supports the theory 
that the strength of the UV upturn in elliptical galaxies is directly tied to that fraction of a 
high Zevoi and high Yzams population evolving along AGBM tracks from the blue end of 
the horizontal branch. Integrations over high Z^voi and high Yzams evolutionary tracks best 
£t the HUT data, and AGBM evolution is more easily produced in a high Z^voi and high 
Yzams population, providing a natural explanation for the correlation between UV upturn 
strength and metallicity found by Burstein et al. (1988). 

Of course, the large uncertainties in our analysis mean that we cannot rule out the Park 
& Lee (1997) hypothesis. The statistics in our htting favor the high Y & Z^voi models, 
but to the eye one can produce reasonable matches to the HUT data with integrations 
over low Y & Zevoi tracks. In the Park & Lee work, they argue that the FUV flux is 
dominated by a population with [Fe/H] < -0.7, but unfortunately it is difficult to discern 
from their paper the exact distributions in M^nv and metallicity which are used in their 
models. Furthermore, although Park & Lee stress that their purpose is not to match the 
actual SEDs observed in ellipticals, it seems that their synthetic spectra are too flat in 
the FUV, at least in comparison to the overplotted lUE data (cf. their Figs. 6 & 10). It 
would be interesting to apply their models to the HUT data, in order to determine if their 
low metallicity models can reproduce the actual SEDs of observed ellipticals, including the 
flux near the Lyman limit, since our analysis indicates that the low Y h Z tracks do not 
reproduce the short wavelength flux needed to reproduce the HUT SEDs. 

In order to answer definitively the debate and disentangle the effects of age and 
metallicity on the HB mass distribution, data with higher S/N for a larger sample of 
galaxies is required. At higher S/N, it should be obvious if diffusion processes are playing a 
large role in the strengths of the FUV absorption lines, since observations of evolved stars 
show lines that deviate strongly from solar abundances or even scaled solar abundances. 
However, we should not restrict our study to the strengths of metallic lines. The hydrogen 
lines are also very sensitive to the effective temperature and surface gravity of the different 
HB evolutionary phases. Inspection of Fig. 1 shows that PAGB spectra have a much 
shallower Lyman series than that seen in AGBM spectra. Although the S/N in the HUT 
data was not sufficient to use the Lyman series as a stellar population probe, accurate 
measurements of the flux below Lyman a (and thus beyond the detection capabilities of 
HST) will yield important information about the relative numbers of PAGB and AGBM 
stars in a given population. Such observations could be made with another flight of HUT 
or a similar successor. Fortunately, future FUV measurements of elliptical galaxies will be 
accompanied by the rapid progress currently seen in the calculation of synthetic stellar 
spectra. The near future should bring the emergence of fully line-blanketed, non-LTE model 
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atmospheres and synthetic spectra, and with these models, another source of uncertainty 
will be removed from our understanding of the FUV data. 

Once we know the chemical composition of the FUV sources in elliptical galaxies, the 
resolution of the debate will have impact beyond the study of stellar evolution. Elliptical 
galaxies show great promise as standard candles and tracers of cosmic evolution. The 
understanding of the FUV flux from galaxies will also enable us to decode the information 
stored in the longer UV wavelengths (Dorman & O’Connell 1996 ). The different 
composition hypotheses each argue for different ages of the Universe and different scenarios 
for the construction of the galaxies. Under the Park & Lee (1997) low-metallicity scenario, 
the more massive galaxies formed first, in contrast to the “bottom-up” cosmologies wherein 
the giant elliptical galaxies formed from the merger of smaller galaxies. In this way, the 
metallicity debate has consequences that extend to important issues in cosmology and 
galaxy evolution. 
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Fig. 1.— The three classes of evolution for stars on the HB produce very different FUV 
spectra when one integrates the FUV luminosity over these tracks, assuming that stars 
enter each track at the rate of one star per year. The gap between the AGB evolution 
{Menv = 0.231 Mq; dashed; DR093) and the PAGB evolution {Mcore= 0.569 Mq; dashed; 
Vassiliadis & Wood 1994) is due to the uncertainties in the TP phase. All of the evolutionary 
paths assume solar abundances, and the luminosities are integrated over the synthetic spectra 
grids of Brown, Ferguson, & Davidsen (1996) and Kurucz (1992), again assuming solar 
abundances. Note that the integrated luminosity from the AGBM track {Menv = 0.006 Mq; 
thin solid; DR093) is much higher than that from the PEAGB track {Menv = 0.051 Mq; 
dotted; DR093), which in turn is higher than that from the PAGB track (dashed). The 
heavy solid line denotes the ZAHB. 









Fig. 2.— The optical image (top panel) and the FUV image (bottom panel) of M 60 
and NGC 4647 demonstrate that there is no evidence for tidal interaction between the two 
galaxies, that much of the concentrated FUV light from M 60 fell within the HUT slit (white 
rectangle), and that light from NGC 4647 is unlikely to contaminate the HUT data for M 60. 
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Fig. 3.— Fig. 3- The best fit evolutionary model (solid) from integrations over EHB tracks 
is plotted over the HUT data (solid; histogram; binned by 2.5 A) for each of the Astro-2 
observations. The evolutionary models integrate over synthetic spectra with Zatm = 0.1 Zq, 
which best reproduces the absorption features in the Astro-2 data, especially in the higher 
S/N data of M 60. Statistical errors in the data are shown dotted. 
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Fig. 4.— The best fit evolutionary model (dashed) from each of the DR093 abundance 
subsets (Table 2) shows that EHB stars with high Zf.voi and high Yzams fit the M 60 data 
(solid histogram) better than the EHB stars with low Z^voi and low Yzams- The abundances 
and envelope mass driving the evolution for each of the evolutionary tracks is given in each 
panel, and Zatm = 0.1 Zq. Part of the reason the models with supersolar evolutionary 
abundances (right column) £t better than those with subsolar abundances (left column) is 
the flux deficiency shortward of 970 A (AF). 
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Fig. 5.— The best fit evolutionary model (solid) from two-component models integrating 
over EHB and PAGE tracks is plotted over the HUT data (solid histogram; binned by 
2.5 A) for each of the Astro-2 observations. As in Fig. 3, the evolutionary models have 
Zatm = 0.1 Z©. Instead of statistical errors, the PAGE contribntion to each two-component 
model is shown dashed. 
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Fig. 6.— The same as in Fig. 4, except that the model (dashed) now has two components: 
EHB stars and PAGE stars. The PAGE component comes from the Zevoi = Zq set of 
Vassiliadis and Wood (1994), with Mcore = 0.569 M©. The legend on each of the eight plots 
reflects the abundances which drive the evolution of the DR093 EHE tracks. The models 
with subsolar evolutionary abundances (left column) again show signihcant flux dehciencies 
shortward of 970 A (AF). 
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Fig. 7.— The (m<iooo~''^i 55 o) colors in the integrated synthetic spectra change considerably 
as a function of envelope mass on the HB (right panel). PAGE stars evolving from the red 
end of the HB and EHB stars evolving from the extreme blue end of the HB both show hotter 
(^^<1000 — ml550) colors, while stars evolving between these extremes appear cooler. Thus, 
in the color-color diagram (left panel), galaxies with strong UV-upturns (NGC 1399, M 60, 
M 89) may be dominated by EHB stars evolving from narrow ranges of envelope mass on the 
blue end of the HB, while galaxies with weaker UV-upturns (M 49, NGG 3115, NGG 3379) 
may have signihcant FUV contributions from PAGE stars. Error bars show the la statistical 
uncertainty, and the dashed lines show the shift in (rriisso — V) color that would occur if we 
used the Burstein et al. (1988) measurements. 
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Fig. 8.— The HUT data for M 60 (histogram) are plotted with the best fit EHB + PAGE 
model (solid) for M 60, given in Table 3. The EHB + PAGE model integrates over synthetic 
spectra with Zatm = Zq, 0.1 Zq, and 0.01 Zq. The absorption features £t best in the 
center panel, with Zatm = 0.1 Zq. Small differences in the data shown in each panel are due 
to variations in the airglow subtraction, since airglow profiles were fit simultaneously with 
the models. 





Metals Present 

Fig. 9.— The spectral indices for synthetic spectra with Tgff = 25,000 K and log g = 5.0 
are plotted for spectra at solar abundance (labeled All) or at solar abundance with one 
element missing (labeled). The line blanketing in the FUV prevents the dehnition of a 
“clean” spectral index that is affected by only one element, but these indices show the most 
sensitivity to the intended elements. The indices are dehned in Table 2. 
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Fig. 10.— The “UV-strong” spectrum (top panel) is the variance-weighted average of M 60 
and M 89, and the “UV-weak” spectrum (bottom panel) is the variance-weighted average 
of M 49, NGC 3115, and NGC 3379. The line and continunm regions used in the spectral 
index measurements are indicated. 
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Fig. 11.— The spectral indices from Table 4, plotted for EHB+PAGB models and the 
Astro-2 data. The UV-weak galaxies (M 49, NGC 3115, and NGG 3379) were averaged to 
form a composite UV-weak spectrum, and the UV-strong galaxies (M 60 and M 89) were 
averaged to form a composite UV-strong spectrum. In order to account for spectral index 
variations due to evolutionary assumptions, we have plotted in each panel a given spectral 
index vs. Zatm for eight EHB-I-PAGB models (one from each of the eight subsets in 
Table 2). Specifically, each EHB-I-PAGB model is taken from the best-fit models to M 60 
that are shown in Fig. 6. Note that the data are consistent with low atmospheric abundances 
(Zatm), with the UV-strong galaxies tending toward stronger indices. 
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Table 1: HUT Astro-2 Observations 



M 49 

M 60 

M 87 

M 89 

NGC 3115 

NGC 3379 

Exp. (sec) 

1346 

5224 

950 

1682 

1634 

3074 

D (Mpc) 

21.9® 

21.9® 

21.9® 

21.9® 

lO.O^^ 

18.2® 

V (mag)“ 

11.19 

11.46 

12.30 

11.60 

10.54 

11.24 

H I (10^0 cm2)d 

1.64 

2.41 

2.51 

2.50 

4.44 

2.86 

E{B - V) (mag)® 

0.000 

0.010 

0.023 

0.035 

0.025 

0.013 

(mi550 - V) (mag)f 

3.60 ±0.03 

2.43 ±0.01 

1.61 ± 0.02 

2.69 ± 0.02 

4.09 ± 0.02 

4.07 ± 0.02 

(m<iooo - mi55o) (mag)s 

-0.21 ±0.05 

0.19 ±0.02 

-0.29 ±0.05 

0.07 ± 0.05 

-0.00 ±0.05 

-0.31 ±0.04 


“Ferguson & Sandage 1990 
^Hanes & Harris 1986 

“Through the HUT slit, from archival HST data. 

'^Einstein On-Line Service, Smithsonian Astrophysical Observatory 
“Burstein & Heiles 1984 

■^Through the HUT slit, from HUT and archival HST data. 
^Through the HUT slit. 


Table 2: Evolutionary Models 


Set 

[Fe/H] 

|0/Fe] 

^evol 

Yzams 

Yhb 

Mcore 

A 

-2.26 

0.50 

0.0001 

0.235 

0.245 

0.495 

B 

-1.48 

0.63 

0.0006 

0.236 

0.247 

0.485 

C 

-0.47 

0.23 

0.0060 

0.238 

0.257 

0.475 

D 

0.00 

0.00 

0.0169 

0.270 

0.288 

0.469 

E 

0.39 

0.00 

0.0400 

0.270 

0.292 

0.464 

F 

0.58 

0.00 

0.0600 

0.270 

0.289 

0.458 

G 

0.43 

0.00 

0.0400 

0.340 

0.356 

0.454 

H 

0.71 

0.00 

0.0600 

0.450 

0.459 

0.434 



Table 3: Stellar Population Modelling- EHB Fitting 



M 49 

M 60 

M 87 

M 89 

NGC 3115 

NGC 3379 

NGC 1399 

SEFtot (stars yr l)®- 

0.94 

0.76 

0.36 

0.72 

0.39 

0.65 

0.84 

Best group 

H 

G 

H 

G 

G 

H 

H 

Menv 

0.086 

0.036 

0.076 

0.036 

0.036 

0.076 

0.016 

x^li' 

1.2 

2.9 

1.7 

1.7 

1.3 

1.3 

3.0 

SEFehb (stars yr-^)^ 
Ax^ for alternative 

2.1x10-2 

5.7x10-2 

5.1x10-2 

4.3x10-2 

6.2x10-3 

9.7x10-3 

8.1x10-2 

abundance groups? 

A 

39 

371 

64 

47 

43 

47 

211 

B 

30 

370 

51 

50 

40 

37 

142 

C 

20 

218 

40 

25 

26 

35 

97 

D 

21 

181 

40 

27 

29 

35 

81 

E 

17 

34 

34 

3 

9 

31 

71 

F 

24 

176 

6 

18 

21 

18 

65 

G 

4 

0 

11 

0 

0 

15 

40 

H 

0 

88 

0 

13 

10 

0 

0 


“Through the HUT slit, based on the fuel consumption theorem and U-magnitude from Table 1. 

^Through the HUT slit, based upon fitting the HUT FUV data. 

“By statistical comparison to the best-fit EHB model for each galaxy, v is 272 for NGC 1399, and 326 for the 
Astro-2 galaxies. 



Table 4: Stellar Population Modelling- EHB + PAGE Fitting = 0.569 Mq) 



M 49 

M 60 

M 87 

M 89 

NGC 3115 

NGC 3379 

NGC 1399 

SEFtot (stars yr 

0.94 

0.76 

0.36 

0.72 

0.39 

0.65 

0.84 

Best Zsvol group 

G&H 

E 

H 

D&F 

B&F 

C&H 

H 

hdenv 

0.026 

0.026 

0.056 

0.017 

0.017 

0.116 

0.016 


1.2 

2.6 

1.6 

1.7 

1.3 

1.2 

3.0 

4f7 limits^ 

0.007-0.116 

0.021-0.046 

0.016-0.106 

0.012-0.046 

0.012-0.166 

0.007-0.166 

0.011-0.026 

SEFehb (stars yr“i)'' 

2.1x10-2 

5.6x10-2 

3.9x10-2 

4.1x10-2 

4.7x10-3 

4.8x10-3 

6.8x10-2 

SEFpAGB (stars yr-l)“ 

for alternative 
abundance groups^ 

0.38 

0.76 

0.36 

0.72 

0.18 

0.41 

0.84 

A 

5 

131 

47 

13 

3 

1 

211 

B 

3 

118 

43 

15 

0 

1 

142 

C 

1 

37 

30 

4 

2 

0 

97 

D 

2 

16 

32 

0 

3 

1 

81 

E 

2 

0 

24 

2 

2 

3 

71 

F 

2 

21 

7 

0 

0 

2 

65 

G 

0 

3 

6 

8 

5 

4 

40 

H 

0 

144 

0 

27 

7 

0 

0 


“Through the HUT slit, based on the fuel consumption theorem and U-magnitude from Table 1. 

^In any Z^voi group. 

“Through the HUT slit, based upon fitting the HUT FUV data. 

'^By statistical comparison to the best-fit EHB-f-PAGB model for each galaxy, v is 271 for NGC 1399, and 
325 for the Astro-2 galaxies. 
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Table 5: Stellar Population Modelling- EHB + PAGE Fitting = 0.597 Mq) 



M 49 

M 60 

M 87 

M 89 

NGC 3115 

NGC 3379 

NGC 1399 

SEFtot (stars yr 

0.94 

0.76 

0.36 

0.72 

0.39 

0.65 

0.84 

Best Ze„oi group 

G 

G 

H 

E 

F 

H 

H 

Menv (-^q) 

0.026 

0.046 

0.056 

0.021 

0.017 

0.076 

0.021 


1.2 

2.7 

1.6 

1.7 

1.3 

1.2 

3.0 

Menu 4cr limits^ 

0.011-0.106 

0.046-0.046 

0.016-0.106 

0.012-0.116 

0.017-0.166 

0.011-0.86 

0.011-0.026 

SEFehb (stars yr”!)'' 

2.1x10-2 

5.3x10-2 

5.1x10-2 

4.1x10-2 

4.9x10-3 

7.5x10-3 

9.7x10-2 

SEFpAGB (stars yr“l)'' 

0.018 

0.76 

0.36 

0.71 

0.33 

0.38 

0.84 

for alternative 








abundance groups^ 








A 

30 

233 

55 

22 

7 

20 

211 

B 

26 

217 

47 

23 

5 

17 

142 

C 

10 

130 

35 

13 

3 

5 

97 

D 

15 

72 

36 

3 

4 

40 

81 

E 

3 

39 

28 

0 

1 

12 

71 

F 

20 

107 

4 

9 

0 

24 

65 

G 

0 

0 

8 

3 

3 

1 

40 

H 

3 

123 

0 

16 

5 

0 

0 


“Through the HUT slit, based on the fuel consumption theorem and U-magnitude from Table 1. 

^In any Z^voi group. 

“Through the HUT slit, based upon fitting the HUT FUV data. 

'^By statistical comparison to the best-fit EHB-I-PAGB model for each galaxy, v is 271 for NGC 1399, and 
325 for the Astro-2 galaxies. 


Table 6: Spectral Indices 


Ion 

Line region (A) 

Continuum regions (A) 

Plot symboP 

cm 

1171-1179 

1161-1171, 1179-1186 

□ 

CII 

1332-1338 

1326-1332, 1338-1347 

A 

N III 

987-998 

981-987, 998-1007 

□ 

N II 

1080-1090 

1074-1080, 1090-1094 

A 

Si III 

1105-1114 

1100-1105, 1114-1120 

□ 

Sill 

1257-1268 

1251-1257, 1268-1275 

■k 

Si IV 

1390-1406 

1380-1390, 1406-1411 

A 

Fe III 

1508-1519 

1502-1508, 1519-1523 

□ 

Fe III, Fe II 

1623-1637 

1615-1623, 1643-1648 

k 


In Figs. 9 and 11. 
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Table 7: Absorption Lines 


Feature 

Ao 

M 60 
EW(A)^ 

Significance 

NGC 3379 

EW (A)^ Significance 

Model 

EW(A)° 

Lyd 

949.8 

2.2 ±0.5 

7.1x10-^ 

2.5 ±0.7 

4.1x10-1 

3.8 

Ly7 

972.9 

3.2 ±0.6 

3.9x10-® 

2.4 ±0.5 

1.1x10-1 

6.8 

cm 

977.2 

1.4 ±0.6 

7.6x10-1 

2.8 ±0.5 

7.2x10-1 

1.6 

N III 

991.7 

< 1.5 

- 

0.8 ± 1.1 

9.2x10-1 

0.8 

Ly/5 

1026.0 

10.6 ±0.4 

< 10-1° 

8.9 ±0.4 

2.0x10-° 

12.7 

CII 

1036.4 

1.1 ±0.4 

2.8x10-1 

0.3 ±0.6 

1.0x10° 

1.5 

N II 

1085.3 

1.5 ±0.4 

3.2x10-1 

0.9 ±0.7 

8.0x10-1 

0.8 

cm 

1175.8 

1.3 ±0.4 

3.0x10-® 

2.2 ±0.7 

8.2x10-2 

0.8 

Si II 

1263.0 

< 1.0 

- 

< 2.7 

- 

1.4 

CII 

1335.2 

0.6 ±0.3 

5.5x10-1 

< 2.5 

- 

0.6 

Si IV 

1394.1 

1.1 ± 1.1 

8.7x10-2 

0.3 ±0.4 

0.99x10-1 

0.3 

Si IV 

1402.5 

1.2 ±0.6 

6.5x10-2 

< 2.0 

- 

0.1 

CIV 

1550.2 

< 1.4 

- 

< 2.7 

- 

0.2 


“Uncertainties are from the Icr errors, and upper limits are for the 2 ct confidence level. 
^Best fit EHB + PAGB model for M 60, with Zatm = 0-1 Zq. 



